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Abstract. We present far-infrared and millimetric observations of a sample of 3C objects. Millimetre data were 
taken at 1.25 mm with the IRAM 30m antenna feeding the MPIfR bolometer array. Mid-infrared (MIR) and 
far-infrared (FIR) photometry were carried out with the ISOCAM and ISOPHOT cameras on the ISO Satellite. 
Additional FIR IRAS observations are also included. We present the entire Spectral Energy Distributions (SEDs) 
from the UV to radio and discuss the emitting mechanisms. Two composite spectra, one for the radio galaxies 
and one for the radio quasars, are built in the object rest frame. While the SEDs of the two classes differ, they are 
indistinguishable in the MIR and FIR range where they probably arise as thermal emission from a dusty torus 
and a larger-scale (cooler) dust distribution in the host galaxy. 
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1. Introduction 

While classical papers on colours and luminosities of 
radio galaxies have explicitly ignored the presence of 
dust on the evidence of the strong Ly-a emission, blue 
colours and small scatter in the K-band Hubble diagram 
( [Lilly fc Lonagir 1984| ; [Punlop et al. 1989| ; [Lilly 19891 ), 



Cimatti fc di Serego Alighicri 1995 



Manzini & di Serego Alighieri 1996; 
Best fc Longair 199S| ; [Vernet et al. 2001D 



more recent work adopts a more cautious approach 
and suggests that some of the observed properties can 
be interpreted as evidence for the presence of dust 
distributed on scales of tens of kpc. For instance, 
the detection of significant linear polarization in the 
UV/blue aligned light is identified with scattered light 
from a hidden AGN due to externally illuminated dust 
which acts as a very efficient reflector of UV light 



Images taken with HST show evidence for dust lanes in a 
large fraction of SCR galaxies with redshift 0.1 < z < 0.5 
and significant extinction is found in the central few 
kpc in some radio galaxies . Dust disc orientation 
is close to being perpendicular to the radio source 
axis and obscuration from dust may dominate the 
appearance of intermediate redshift sources, influence 
the apparent morphology of the optical galaxy and 
determine the alignment effect seen in powerful ra- 



dio galaxies ( 


Baum et al. 1995; 


de Koff et al. 1996 


McCarthy et al. 1997 


; dc Koff et al. 2000). 



([di Serego Ahghieri et al. 1989 


; Bcarrott et al. 199C 


Cimatti et al. 1993 


; di Serego Alighieri et al. 1996 
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Furthermore, comparison of line strengths arising from 
the same ions and involving a common ground state 
with calculated line ratios provide a sensitive measure 
for reddening. Attempts were made to measure in a 
number of high-z radio galaxies and, in most of the 
cases, a significant amount of reddening (with Ay > 0.3) 
is found. 

On purely theoretical grounds, one expects the produc- 
tion of significant amounts of dust at early epochs when 
galaxies were undergoing vigorous star formation. Thus 
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there is reason to believe that much of the UV luminos- 
ity of high redshift radio sources is reprocessed by dust. 
The question of whether there is indeed a large amount of 
dust associated with radio galaxies is important not only 
from the point of view of galaxy formation and evolution, 
but it can help us understand the apparent differences be- 
tween radio galaxies and radio quasars and as a test for 
the Unification Scheme for AGN (van Bemmel, Barthel, 
de Graauw, 2000). 

Direct emission from dust can be detected at FIR 
and sub- mm wavelengths where the thermal re-radiation 
from dust grains is expected. However, on the basis of 
IRAS data alone for quasars and radio galaxies , it is 
not straightforward to constrain the emission mechanism 
with great confidence and, in particular, the relative con- 
tribution of thermal and non-thermal components at FIR 
wavelengths can be addressed only in sources which are 
strongly lobe-dominated and therefo re, supposedly, free of 
any observable beam ed radiation ( Heckman et al. 1992 
H eckman et al. 1994| ; 



Hes et al. 1995 



Hockstra et al. 1997; van Bemmel ct al. 2000). In these 



objects the radio axis is further away from our line of 
sight, and the beaming of non-thermal emission therefore 
less efficient. Thus if dust is present, one might see it in 
these objects. 



The detection of HI absorption (Uson et al. 1991 



Carilli et al. 1998) and CO emission has confirmed the 
existence of dense concentrations of both atomic and 
molecular gas in the nuclei of some radio galaxies and 
quasars and indicated that rich supplies of molecular gas 
may be ubiquitous in powerful radio objects detected 
by IRAS ([Mirabel et al. 1989t [Mazzarella et al. 1993 ; 



Evans et al. 1999| ; [Evans ct al. 200lT 



Detection of sub- mm continuum (with JCMT at 850 /im 
and with IRAM 30m at 1.25 mm) in radio galax- 
ies provides a picture in which sources with typical 
redshift z < 1 are, on average, not detected, while 
those at higher redshifts (z > 2.4) present strong sub- 
mm fluxes, suggesting the presence of spatially dis- 
tributed dust with approximately 10^ Mq in these objects 



(Chini & Kriigel 1994; 


Dimlop et al. 1994 




[vison 1995 


Best et al. 199^; 


Ivison et al. 1998 




van der Werf 1998 


Archibald et al. 2001 


). On the one 


hand, this could sim- 



ply reflect a K-correction effect: the steep (S oc 



,3-4 



slope of the dust-emitted spectrum producing an appar- 
ent constancy or increase of observed flux, at a flxed ob- 
served frequency, with redshift at constant luminosity. 
Alternatively, one could envisage an evolutionary effect 
with larger dust masses and luminosities at higher red- 
shifts. In the latter case, the large FIR luminosities could 
be produced by large star-formation rates in extremely 
gas-rich galaxies, but it can be also due to a selection effect 
which favours the d etection of those h i gh-z objects wit h 
greater radio power ( van der Werf 1998 ; Best et al. 1998[) . 

At low-redshifts, the 3CR sample was the subject of 
extensive investigation at FIR wavelengths but the ques- 
tion of the presence of dust in these objects could not be 
settled on the basis of IRAS data alone. 



Heckman et al. (1992,1994) claim that, if long wavelength 
FIR emission is due to thermal re-radiation by circum- 
nuclear dust, quasars and radio galaxies are expected 
to show similar outputs of such radiation, because of the 
optical thinness of the thermal emission. But they show 
that 3C quasars are more powerful FIR (60 /im) sources 
than 3C radio galaxies and this result can be used either 
to disprove the uniflcation scheme of quasars and radio 
galaxies proposed by Barthel (1989), or to state that the 
FIR emission from quasars cannot be due to thermal rera- 
diation from dust unless the emitting region is very opti- 
cally thick, resulting i n some degree of radiated anisotropy 
( |Pier fc Krohk 199^ [Granato fc Danese 1996] )). 

This issue was the subject of subsequent investi- 
gations. By computing the relative contributions from 
a relativistically-beamed core and isotropic emission at 
60 /im, a signiflcant beamed component from the nuclei 
of lobe-dominated quasars and radio galaxies in addi- 
tion to the isotropic thermal dust component was found 
( |Hes et al. 1995[ [Hoekstra et al. 1997| ). The infrared out- 
put of radio galaxies and even some quasars should also 
be affected by contributions from star-formation episodes 
(van Bemmel et al. 1998). 

Deep integrations performed with the ISOPHOT camera 
on the ISO Satellite on a small sample of radio galax- 
ies and quasars conflrm the previously reported FIR 
excess in quasars (van Bemmel et al. 2000), with this 
excess extending up to restframe wavelengths of 130 
/im. In their radio galaxies the relative contribution of 
dust at FIR wavelengths with respect to beamed emis- 
sion was estimated to be 98% and emitted by cool dust. 
Meisenheimer et al. (2001) on the basis of 13 detections 
out of 20 ISOPHOT observations on 10 pairs of 3C radio 
galaxies and quasars do confirm the thermal nature of 
the infrared emission, but do not conflrm the presence of 
an infrared excess in quasars. 

We present in this paper ISO FIR and groundbased 
millimetric observations of a sample of 3C sources (see 
Sect. 2), containing both quasars and radio galaxies . 
We discuss the origins of the FIR- mm emission and the 
differences found between the two classes. Data are pre- 
sented in Sect. 3 and the results are discussed in Sect. 
4. 



2. The sample selection 

The sample discussed in the present work was selected 
from a larger sample which was proposed (proposal ref- 
erence MRC-3CR) for observation with the ISOPHOT 
(Lemke et al. 1996) and, for a subset of sources, with 
the ISOCAM (Cesarsky et al. 1996) cameras on the ISO 
Satellite. The 3C sources were selected solely on the basis 
of their visibility with the ISO satellite in regions of low 
Galactic cirrus emission. No objects with z < 0.3 were 
included and the sample contained 50 sources classified 
as quasars and 85 as galaxies. The main purpose of the 
ISO proposal was an extensive study of the FIR emission 
mechanism in 3C sources and a comparative study of ra- 
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dio galaxies and quasars . These data, supplemented 
with millimetric observations, were to be used to disen- 
tangle radio synchrotron emission from the thermal dust 
emission. The assembled sample spans a large range of 
radio properties and cosmic look-back times to avoid the 
well known biases in high frequency selected samples that 
arise from relativistic beaming effects, while maintaining a 
statistically useful mix of radio galaxies and quasars. This 
strategy demands complete samples selected at low fre- 
quencies without any biases in spectral index, radio mor- 
phology or angular size. The SCR sample is well suited 
for this problem within the redshift range that it spans 
(z < 2) and is still the only completely optically identified 
sample of low-frequency selected objects. This sampling 
frequency (178 MHz) does guarantee selection due to un- 
beamed, optically thin and thus isotropic, emission from 
the radio lobes and so avoids introducing any orientation 
bias. 

After launch, the sensitivity of ISOPHOT proved to be 
inadequate to complete observations effectively in the time 
available. In the present work we discuss the observations 
taken with the IRAM 30m antenna, feeding the MPIfR 
bolometer array at 1.25 mm, of 27 of the 3C sources, se- 
lected purely on the basis of telescope visibility, from the 
original ISO proposal. We also present the ISOPHOT ob- 
servations in the wavelength range 5-200 /im of 15 of them 
and ISOCAM observations of 10. Even though the origi- 
nal ISOPHOT observations were not completed, it is still 
possible to address some of the purposes of the proposal 
in a self-consistent way. 

3. Observations 

3.1. 1.25 mm data 

The 1.25 mm data presented here were taken with the 
MPIfR 19-channel bolometer ( |Kreysa et al. 1998D at the 
focus of the IRAM 30m antenna (Pico Vcleta, Spain) 
during March 1996 and March 1997. The filter set com- 
bined with the atmospheric transmission produces an ef- 
fective wavelength around 1.25 mm; the beam size is 11" 
(FWHM) and the chop throw was set at 32" with a 
chopping frequency of 2 Hz. The average sensitivity for 
each channel, limited principally by atmospheric noise 
and measured before any sky-noise subtraction was 60 
mJy/Vffz. The effect of the sky noise on flux measure- 
ments could, however, be substantially reduced by exploit- 
ing the correlation between signals from the different chan- 
nels using the standard three beam (beam-switching -|- 
nodding) technique. The average rms value was 1 mJy for 
a typical integration time on-source of 2000 s. 

Atmospheric transmission was monitored by making 
frequent skydips from which the derived zenith opacities 
were 0.09-0.3. Calibration was performed using Uranus as 
primary calibrator and Mars and Quasars from the IRAM 
pointing list as secondary sources. The different measure- 
ments vary by less than 5 % for both planets. If we include 
the uncertainty in the planet temperature, we estimate an 



average flux calibration uncertainty of 10 %. Pointing was 
checked each hour and the average accuracy achieved was 
better than 3". 

The data were reduced assuming that the target 
sources are unresolved, i.e., having an extent at mm wave- 
lengths smaller than the size of the central channel. The 
other 18 channels (excluding one which suffered a large 
electronic loss) were then exploited to derive a low-noise 
sky estimate. The weighted average value of the sky, com- 
puted using these outer 17 channels, was subtracted from 
the signal in the central channel. Note that this procedure 
eliminates only that part of the sky fluctuation with corre- 
lation length smaller than the chop throw (32"), i.e., from 
fluid motions at short wavelengths. However, the domi- 
nant part of the atmospheric noise is produced by motions 
of convective shells on large scales (at long wavelengths) 
(Church, 1995), while high frequency (5-20 Hz) fluctua- 
tions do not contribute significantly to the noise and are 
averaged over 0.25 s per phase, while the wobbling of the 
secondary smooths out the low frequency noise. According 
to Andreani et al. (1990), the correlation length for con- 
vective shells at these wavelengths is of the order of sev- 
eral tens of centimeters, i.e., only fluctuations generated 
at an altitude greater than 2000m above the telescope 
survive the double-switching subtraction and contribute 
to the noise. At altitudes greater than 5000m above sea 
level, however, the residual water vapour is very low and 
contributs little to the noise. 

1.25 mm fluxes with the associated la uncertainties are 
given in table 1. The first uncertainty is statistical while 
the second one derives from the calibration uncertainty 
related to the planetary measurements and their tem- 
perature uncertainty. Our 1.25 mm fluxes towards 3C286, 
3C 309.1, 3C295 and 3C325 differ slightly from those re- 
ported by Meisenheimer et al. (2001). They agree, how- 
ever, within the calibration uncertainties. 



3.2. IRAS and radio data 

The IRAS data are taken from co-added survey data 
provided by IPAC and based on the SCANPI (Scan 
Processing and Integration Tool) processor. This proce- 
dure computes the one-dimensional co-addition of all the 
IRAS survey data of the source. The sensitivity is compa- 
rable to that achieved by the FSC (Faint Source Catalog) 
for point sources (see the IPAC manual for details). The 
resulting fluxes are listed in table 2. We have also ap- 
plied the SUPERSCANPI procedure, which is similar to 
SCANPI but allows the inference of an average value of 
the FIR emission of the whole class of object. By dividing 
our sample into two separate classes, radio galaxies and 
quasars , we have applied the SUPERSCANPI procedure 
to estimate the weighted-average fluxes at 12, 25, 60 and 
100 /im separately for each. These values are listed in ta- 
ble 2. The same procedure, but applied to a much larger 
set of objects, was used by Heckman et al. (1992). Within 
our larger uncertainties our results for radio galaxies and 
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quasars are in agreement. For 3C 268.3 IRAS data were 
taken from Hoekstra et al. (1997). Radio data are taken 
from NED (NASA Extragalactic Database). 

3.3. ISO data 

Fifteen objects of the present sample were observed 
with the ISOPHOT camera. For two objects, 3C 268.4 
and 3C 280, data were collected for our programme 
(RFOSBURY MRC-3CR2) on July 3'''* 1996, at 60 and 
90 /^m. Data for 3C295, 3C 309.1, 3C325 and additional 
observations of 3C280 were taken from 5 to 170 /xm, 
those of 3C286, 3C287 from 5 to 100 /xm by Chini and 
have been independently published by Meisenheimer et 
al. (2001). These data together with those of 3C 268.4 
(at 90 and 170 ^im) 3C313 (90Atm), 3C 288.1 (25-170 /im), 
3C352 and 3C356 (90 and 160 /im) were taken from 
the archive and reduced by us w ith the Phot Intera ctive 
Analysis tool (PIA) version 8.0 ( [Gabriel et al. 1998| ). AU 
data were taken in chopping mode except those for 3C 352 
and 3C 356 which were mapped in raster mode. These 
latter data were reduced as described in van Bemmel et 
al. (2000). Data for 3C46, 3C 268.3, 3C295, 3C337 and 
3C 343.1 are taken from Fanti et al. (2000) and that paper, 
in which there is an extensive discussion, should be con- 
sulted for any details concerning the observing procedure 
and the data reduction. 

Our analysis of the ISO data of the sources 3C286, 
3C287, 30 295, 30 309.1, 30325 agree within the error 
bars with that by Meisenheimer et al. (2001) but clear 
detections towards 30 280 were obained because of the 
different observational setup with much longer exposure 
time. 

Our reduction procedure is briefly summarized as fol- 
lows: 

— Ramps (V/s) of the Edited Raw Data (ERD) are lin- 
earized and cosmic hits are removed using the two- 
threshold deglitching procedure. A linear fit is then 
applied to the ramps providing the Signal per Ramp 
Data (SRD). 

~ All the data corresponding to one chopper position are 
then averaged and extrapolated in time to a value of 
constant detector response {Pattern Analysis). 

— To correct the remaining cosmic hits, another deglitch- 
ing procedure is applied. All the ramp signals per 
raster point are averaged and data are then processed 
to Signal per Ohopper Plateau (SOP) level. 

— Data are then corrected for reset interval and dark 
current is subtracted. A signal linearization which cor- 
rects for a changing response during the measurements 
is then applied. 

— The same procedure is applied to the Fine Oalibration 
Source (FSO) data to provide the internal calibration 
data. 

— Source signal is estimated by subtracting the back- 
ground. This latter is determined by averaging the sig- 



nals corresponding to off-source positions of the chop- 
per. 

The resulting fluxes are listed in table 3. Note that for 
30 352 and 30 356, although our analysis shows a signal- 
to-noise ratio larger than 7 at 160 /im and larger than 5 
at 90 /^m, respectively, a good estimate of the flux lev- 
els cannot be achieved because of excessively noisy FSO 
data, which prevent us from performing a good internal 
calibration. 



3.4. ISOCAM data 

ISOOAM observations were taken towards 30 288.1, 
30295, 30 305, 30 309.1, 30 343.0, 30 343.1, 30345, 
30352 and 30356. These ISOOAM observations together 
with other radio galaxies and data reductions will be 
published elsewhere (Siebenmorgen & Freudling, 2001, in 
preparation). A complete ISOOAM catalogue was built by 
these authors. They have reduced staring and raster ob- 
servations of 30 sources in the ISOOAM archive using a 
homogeneous procedure. The procedure was optimized for 
faint sources, and particular effort was taken to effectively 
remove 'glitches' in the data. Aperture Photometry was 
carried out on all sources using two independent proce- 
dures for background estimates and weighting. We address 
the reader to that paper for all further details. 

3.5. UV, optical, near-IR and sub- mm data 

Optical and UV and near IR photometry are taken from 
Best et al. (1998) deVries et al. (1998) and and deKoff et 
al. (1996). Two objects, 30322 and 30 356 were observed 
with SOUBA at 450 and 850 /im by Archibald et al. (2001). 

4. The Spectral Energy Distribution of Radio 
Galaxies 

Fig.s 1-3 show the spectral energy distribution from UV 
to radio wavelengths of the sample sources. 17 out of 27 
objects were detected in the MIR and/or in the FIR. We 
consider as reliable detections all the MID and FIR fluxes 
with > 5, and as marginal detections those with lower 
§ ratio (3 < ^ < 5). 

4.1. Radio spectrum 

We first consider the radio region of the SED. There is no 
way to constrain uniquely the non-thermal radio spectrum 
unless high-frequency {v > 15 OHz) radio data are avail- 
able. A break frequency, depending on the magnetic field, 
for a synchrotron spectrum is expected at high frequencies 
because of electrons cooling. The amount of steepening 
cannot be estimated with the present data since most of 
the sources were not observed in the frequency range be- 
tween 5 GHz and 240 GHz (1.25 mm) and only one third 
have data at 15 GHz. It is, therefore, not possible to re- 
liably extrapolate the radio data at higher frequency to 
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Fig. 1. Spectral Energy Distributions (SEDs) from radio to UV wavelengths, mm points are shown as open circles, 
ISO data as filled circles, IRAS points at 12, 25, 60 and 100 /^m are shown as asterisks, radio and optical observations 
as open squares. For fluxes with signal-to-noise ratios, < 3 a Sfi upper limit is shown as a down-arrow. Two simple 
approximations of the non-thermal emission spectra are shown through the radio data: the solid curve corresponds 
to equation || while short-dashed lines to the classical power-laws with spectral index of arg = —1.04 for radio 
galaxies and aqs — —0.91 for quasars. The long-dashed curve is the output from a model whose basic prescription is 
the reprocessing by ISM dust of the UV-optical light (see text for details). 



disentangle any thermal contribution from the 240 GHz 
flux and infer or reject the presence of another component 
contributing to the mm emission. For the 17 objects de- 
tected in the FIR, it is possible to combine the FIR data 
with the mm point and try to estimate the two likely 
contributors — thermal and non-thermal. 

In what follows, we make two different assumptions. 
First we assume that the sources are all sufficiently old 
that the turn-over frequency is low. In this case, most of 
the electrons have lost their energy since there is no mech- 



anism to continuously produce and/or accelerate them. 
We only fit points at frequencies above the turn-over (as- 
sumed to be at 178 MHz) with the usual power-law of 
synchrotron emission F oc v". We take as a the average 
values for 3C radio gal axies and quasars , i.e . arg = 1-04 



■*qso 



0.91 found by ( Heckman et al. 1992 ) in the fre 



quency range 1.4-15 GHz. We then assume that the syn- 
chrotron spectrum between a few GHz and hundreds of 
GHz maintains the same slope and we extrapolate it to 
mm wavelengths. Any difference between the predicted 
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log i/(Hz) 

Fig. 2. same as figure 1 



and observed values at 240 GHz is attributed to an addi- 
tional spectral component. 

This simple approximation provides a poor match to 
the whole radio spectrum, low and high frequency radio 
data lie above and/or below the expected power-law. In 
particular, the mm emission for all radio galaxies and 
the three quasars 3C268.4, 3C 280.1, 5C 5^5. i lies well be- 
low this extrapolation (see also van Bemmel & Bertoldi, 
2001). It must be stressed here that the 1.25mm observa- 
tions only refer to the flux emitted by the source within 
the central 11". If the radio- lobes dominate the radio spec- 
trum and have a larger extent, it is possible that some 
1.25 mm flux is lost in our observations. This could be 
true for the giant radio galaxies 3C277 and for 3C356. 
For all the other objects, which have a more compact mor- 
phology, the further steepening of the spectrum is real. 



At present the contribution from radio-lobes to the mm 
flux it is not observationally settled. For instance, van 
Bemmel and Bertoldi (2001) do not detect any difference 
in radio-millimetre SED of large and small objects and 
suggest that at least in their objects is the core dominat- 
ing the millimetre emission. 

More generally, a self-absorbed synchrotron emission 
spectrum with an electron power-law energy distribution 
N(E) oc E~^, can be parameterised as follows (see e.g. 
Polletta et al. 2000): 

oc {-r^ {1 - exp{-{-r^-'^-)} e-^c (1) 

where a\ and a-i are the spectral indices for the op- 
tical thick and thin cases respectively (ai = 2.5 and 

0,2 = ^('S — l)/2 for a homogeneous source), vi is the fre- 
quency at which the plasma optical depth reaches unity 
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Fig. 3. same as figure 1 
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and Vc is the cut-off energy of the plasma energy distri- 
bution. For frequencies lower than vt and by neglecting 
the high-energy cut-off, Vc - which although depending 
on the assumed electron energy and the source magnetic 
field for relativistic electron and a Galactic magnetic field 
it is higher than a few hundreds GHz - this expression can 
be approximated by a parabola-shaped curve in a log-log 
plane: 

\ogF^ = C+^ (log v-logvtf (2) 

where \ = 2a\ — a^. We then use this latter expression 
to fit all the radio data with v < 240 GHz where z^t is 
taken equal for all sources, 31.6 MHz. The free parameters 
are C, the normalization factor, and A representing the 
bending of the curve whose best-fit value varies between 
-0.2 and -0.55. For those objects with more than 4 points 



at frequencies lower than 240 GHz a minimizing-x^ fit was 
carried out, while for those sources with fewer than 4 radio 
points the fit is not statistically significant and must be 
considered only as indicative. 

These parabola-shaped curves are shown in Fig.s 1- 
3 as a solid line, while the power-law curve is shown as 
short-dashed. 

For 3C 268.3, 286, 293.1, 295, 343, 343.1 a parabola- 
shaped non-thermal spectrum fits well the entire radio 
spectrum up to 240 GHz. Again we stress here that in 
most cases and in particular for these objects the 240 GHz 
point undoubtedly shows that the radio spectrum bends 
at high energy. However, it is difficult to disentangle from 
the mm data any residual emission not related to the syn- 
chrotron one but to a thermal cool dust component. Even 
for those objects with well sampled radio spectrum only 
the combination of mm with FIR data allows the infer- 
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ence of a different emitting meclianisms on a more solid 
basis. 

3C 286, 287, 309.1 and 345 have a very large mm flux. 
The overall spectrum from radio to optical wavelengths 
could be interpreted as due to a dominant non-thermal 
component. As expected, it is not possible to fit the entire 
spectrum with a single component, as shown in Fig.s 1-3. 
It is likely that the electrons population does not have 
a single age for the entire source. Furthermore, variabil- 
ity plays a role at least for 3C345. Even for 3C286 and 
3C 309.1 with FIR detections it is tough to infer any ad- 
ditional, most likely thermal, component from cold dust 
contributing to the mm flux. 

4.2. FIR and optical spectrum 

For wavelengths A < 200/xm, the interpretation of the ob- 
served SED is not straightforward since both stars and 
AGN, each with their associated dust obscuration and 
emission, must be taken into account. To effectively con- 
strain all these components requires high quality data 
between 3 and 300 /im in addition to those in the 
NIR-optical range. Existing observations allow the con- 
struction of useful SED for thirteen sources with FIR 
data: 3C280 (IRAS/ISO), 3C286 (ISO), 3C 288.1 (ISO), 
3C289 (IRAS), 3C295 (ISO), 3C305 (ISO), 3C 309.1 
(ISO/IRAS), 3C322 (IRAS), 3C325 (ISO), 3C 343.1 
(ISO), 3C345 (IRAS/ISO), 3C352 (ISO) and 3C356 
which have reliable measurements in this range. 

Although the number of data available for most of the 
sources is small, we use them to compare the observed 
SED from the optical to 1 mm with the model developed 
by Mazzei and De Zotti (1996) for radio galaxies . They 
have constructed a spectrophotometric population synthe- 
sis model incorporating dust extinction and re-emission 
and a non-thermal central source which successfully re- 
produces the SEDs of high-z radio galaxies. 

No attempt is made to fit this model to the radio 
galaxies data. We do, however, plot it in Fig.s 1-3, arbi- 
trarily normalized at the 60 ^m flux (or upper limit). In 
particular this spectrum was used to reproduce the radio 
galaxy SEDs where the scattered AGN component dom- 
inates the UV rest frame light and the old stellar pop- 
ulation of the host galaxy is the major contribution at 
longer wavelengths. The dust reprocessed starburst light 
dominates the spectrum at long wavelengths (A > 10 /im), 
while the reprocessing of AGN energy is neglible. This is 
consistent with the model in Vernet et al. (2001) for the 
z~ 2.5 radio galaxy, 4C-h48.48. 

Note that the shape of the predicted spectrum agrees quite 
well with the observed behavior of the radio galaxies UV- 
optical-IR SED, in particular for 3C280, 3C 289, 3C295, 
3C305, 3C322, 3C 343.1, 3C352 and 3C356. Quasars, on 
the other hand, have a different UV-optical spectrum, 
brighter by one dex and dominated by the AGN com- 
ponent. 



4.3. The composite spectra 

To better address the relative importance of the different 
components in radio galaxies and quasars , their com- 
posite spectra were built and are shown in Fig. 4. The 
sample was split into two: radio galaxies and quasars. 
3C 345 was not included in the latter class because of 
its extreme variability. In order to compare spectra of 
sources at different distances data were normalized to the 
60 /im flux. The frequency range is divided into bins in 
the object rest-frame and the 60 /tm-normalized data are 
(weighted) averaged in each bin and errors are given by 
the dispersion around the average. When censored (up- 
per limit) data are present, the Kaplan- Meier estimator 
was used to estimate the average values and their errors 



(Feigelson 1990; LaValley et al. 1992). Here we would like 
to point out the limits of this analysis. It could be argued 
that the mid-far-IR spectra of QSOs do not form a very 
homogenous classs and on the basis of the present data 
alone there seem indeed to be a difference between 268.4, 
280.1, 288.1, 343.0, on the one hand, and 286, 287, 309, 
345, on the other hand, we feel however that it is prema- 
ture to argue for a qualitative difference, because of the 
poor sampling of the spectra of the former objects, and 
we continue to consider the QSO class as a whole. 
Inspection of Fig. 4 allows us to infer the following: 

— The two classes of object show a similar radio spectral 
shape up to a frequency of around 300 MHz. 

— For f > 300 MHz the quasars spectrum is elevated by 
a factor larger than 3 up to 240 GHz (1.25 mm), which 
is most likely to be due to a non-thermal, beamed com- 



ponent (van Bemmel & Bcrtoldi 2001) 



The sharp turn-over of the radio spectrum excludes a 
significant contribution of the non-thermal component 
to the FIR fluxes for the radio galaxies and half of 
the quasars studied in this work. 
No evident difference is detected in the FIR part of the 
spectrum from 100 /im to 1 /xm. This result is satisfac- 
torily in agreement with the fact that the emission has 
a thermal origin, very likely linked to the concurrent 
Starburst in the galaxy disc. This is consistent with 



the result of previous studies ( PoUctta et al. 200C ). 
— The optical and UV part of the spectrum differs by 
one dex or more between the two classes. 

In the framework of Unified Scheme of AGN this re- 
sult is not unexpected since the difference between the 
two classes is only an orientation-dependent effect. An in- 
spection of Fig. 4, together with results from independent 
observations, seems to support this picture: (a) the dusty 
torus emits between 10 and 50 /zm and heated by the AGN, 
(b) an additional cool component emitting at 100-140 /im 
and representing the reprocessed stellar energy in the host 
environment is clearly seen in radio galaxies and it may 
be outshone by the stronger non-thermal beamed com- 
ponent in quasars, (c) the difference at UV-optical wave- 
lengths between the two classes is due to the obscuration of 
the AGN in the radio galaxies by the optically thick torus. 



Paola Andreani et al.: Far-infrared/millimetre Emission in 3C sources 



9 



1000 


1 1 1 1 1 


1 1 1 1 1 1 1 1 


100 




F(i/) (mJy) 


1 






0.1 






0.01 






0.001 






1,0001 


- #050 
r 'RG 




10-= 


r 1 


^] 



• 



1 2 

bg V 



1 6 



Fig. 4. The composite spectra of 3C quasars and 3C ra- 
dio galaxies in our sample, quasars points are shown 
as asterisks, while radio galaxies points as filled circles. 
Data are averaged in each wavelength bin in the object's 
rest- frame (see text for details). Error bars are in general 
smaller than the points size except at log v—\2.'i (150 /im) 
and at logi/=15.25 (0.17/im). 



The optical/UV flux which remains visible in these objects 
can be a combination of scattered AGN light and a direct 



view of a young stellar population (Vernet et al. 2001) 



5. A colour-magnitude diagram for radio galaxies 

Before analysing flux ratios we show in Fig. 5 1.25 mm and 
60/j,m fluxes against the source redshift. Both detection 
and upper limits are equally distributed in redshift and 
do not show any clear trend with z. 

Fig. 6 shows the ratio between the 1.25 mm and the 
60 /xm fluxes, ^^^^ , against the source redshift. 60 /im 
was chosen since it is the most common measured wave- 
length amongst our objects. Asterisks refer to quasars and 
filled circles to radio galaxies. The values for those sources 
with upper hmits at both wavelengths (3C46, 3C 280.1, 
3C 293.1, 3C 305.1, 3C323, 3C337) are identified as open 
circles and plotted at the positions of these hmits. Upper 
limits are given for those objects detected at 60 /im 
but not at 1mm (3C356, 3C277, 3C313, 3C322) and 
lower limits for 3C 268.1, 3C 268.3, 3C 287, 3C 292, 3C 343, 
3C 343.1, 3C 288.1 3C 320, detected at 1 mm but not in the 
FIR. Note that there is no redshift-dependence of the de- 
tection rate at either wavelength. The average value of 
for the radio galaxies class is around 0.05 while 



f 60/xn 



for quasars , it is one order of magnitude higher as seen 
directly from the relative SEDs (Fig.s 1 through 4). A bi- 
modal distribution appears in Fig. 6 where objects whose 
SED is dominated by the non-thermal emission lie in the 




Fig. 5. Upper panel reports the 1.25 mm fluxes, lower 
panel those at 60 /im against the source redshift. Down 
arrows correspond to upper limits. 
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Fig. 6. The ratio """ , is shown against the source red- 
shift. Filled circles are radio galaxies , asterisks quasars . 
Open circles correspond to objects not detected at either 
wavelength. Variation with z of the ratio for a thermal 
spectrum, e(A) • B{X,Td), with Trf^50 K and e cx A^-^ 
(where e is the wavelength-dependent dust emissity) is 
represented by a solid line T^^lOO K by the dotted line 
and Td—20 K by the dashed fine. 

upper part of the plot. In those objects it is likely that any 
thermal mm emission is overwhelmed by the synchrotron 
components. 

The lower part of the plot contains those sources for 
which the thermal component starts to dominate even at 
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Fig. 7. Tlic ratio p'-™'" , is plotted against the source 60 
/Km luminosity. Filled circles are radio galaxies , asterisks 
quasars . Open circles correspond to objects detected at 
neither wavelength. 

mm wavelengths. Some of the objects distribute them- 
selves along a line corresponding to a thermal spectrum 
with a dust mean temperature of 50 K at different red- 
shifts. Objects like 3C295, 3C 268.3, 3C320 and 3C 343.1 
fall to the left of the plot because of a lower temperature 
dust component or to a higher non-thermal component, 
while 3C322, which lies in the lower-right corner could 
have a different thermal spectrum with a hotter average 
dust temperature (around 100 K). This is indeed sup- 
ported by the non-detection of this object in the SCUBA 
photometry at 450 and 850 ;um (Archibald et al., 2001). 

To further investigate the ratio between the 1.25 mm 
and the 60 /xm fluxes, pl^™' , and make it distance- 
independent in Fig. 7 p^""" is plotted against the 60 



fim rest-frame luminosity, L 



for all the sources in 



our sample. The luminosity is computed in an Einstein-de 
Sitter Cosmology with = 1 and Hq = 65km/s/Mpc and 
using a K-corrected thermal spectrum with dust emissiv- 
ity index of /3 = 1.5. 

While the conclusions drawn from this plot should be 
regarded as tentative, the shifting of lower limits, towards 
the left along the x-axis and upwards along the y-axis, as 
indicated by the arrows would not change the observed 
trends. In this figure, the radio galaxies are located pref- 
erentially along a sequence of slightly decreasing p^"™ for 
increasing LeoAtm) while half of the quasars show much 
larger ^^mm values. 

The diagram suggests some interesting conclusions. 
Half of the quasars share the common properties of the 
radio galaxies , while the other half show large p^""" ra- 
tios. This can be explained if these latter quasars have 
a dominant non-thermal spectrum, possibly the beamed 



component much stronger than that observed for those ob- 
jects on the lower-left part of the diagramme. Low value 
of the ratio ^^J""* can be due to the presence of cool dust. 

Most of radio galaxies fall in the region where ( ^ 



Leoiim) is consistent with the expected behavior of a ther- 
mal spectrum with temperatures ranging between T=20K 
and T=50K. Four objects, 280, 356, 289 and 322, lying to 
the right in the plot, could be characterized by a higher 
dust temperature. Indeed the ratio ^^J""* for a thermal 
spectrum with temperature increasing from T=20 K to 
T=200 K at z = 0.5 and z = 1.5 (roughly the redshift 
range of the present sample) decreases as the 60 /xm lumi- 
nosity increases. 



6. Conclusions 

The investigations of the entire spectral energy distribu- 
tions of a small sample of 3C sources with mm and FIR 
observations allows us to draw the following conclusions: 

1. Detected fluxes and upper limits at 60 /zm and 1.25 mm 
are equally distributed in redshift and do not show any 
trend with z. This result is true also for all the ISO 
detections. 

2. The average power-law spectrum through the radio 
data, with spectral indices taken from other studies 
on quasars and radio galaxies (Heckman et al., 
1992), shows that for most of the objects the mm point 
lies below the extrapolation of this law to 240 GHz 
(1.25mm). This means that the radio spectrum bends 
at frequencies lower than 240 GHz. In some of the mm- 
detected objects different spectral components could 
contribute to the detected fluxes. The origin of the 
mm emission in these objects is mainly non-thermal 
for flat-spectrum sources. In some sources there could 
be an additional thermal component from cool dust, 
but without sampling the submm region of their spec- 
trum this component still remains elusive. 

3. QSOs have much stronger mm emission. For half 
of them, the entire SEDs from radio to UV wave- 
lengths can be due to the superposition of differ- 
ent non-thermal components arising from either self- 
absorbed synchrotron and/or synchrotron emission 
from a hard electron spectrum. However for the re- 
maining quasars , the ratio ^^J""* is similar to that 
of radio galaxies , implying a common origin of the 
mm-FIR emission. 

4. Composite spectra for radio galaxies and 
quasars have been constructed in the object restframe: 
the main differences are seen in radio - mm and opti- 
cal SEDs, while the FIR SEDs are remarkably similar. 
This points to a common origin of the FIR emission 
which is likely to be dust-reprocessed energy from star 
formation in the host galaxy and/or with a contribu- 
tion from the dusty AGN-torus. Although on the basis 
of these data alone we cannot fully address this issue in 
quasars, the far-IR detections in quasars suggests that 
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there is a thermal as well as a powerful non-thermal 
component in these objects. 
5. The similarity of the FIR SED of the two populations 
is consistent with the predictions of the orientation- 
based unified scheme and suggests the presence of dust 
either encircling the AGN or in the body of the host 
galaxy or both. On the basis of these data alone it is 
not possible to identify the main heating mechanism 
— AGN or star-formation — though the presence of 
cool dust emission, the width of the FIR spectrum and 
modeling suggest that, at least in radio galaxies , 
dust — which could extend to large distances from 
the heating source — is the main FIR emitter. 
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Table 1. 1.25 mm Observations of 3C sources 



name 


redshift 


type 


Flux 

(mJy) 


stat. unccr. 
(mJy) 


cal. uncer. 
(mJy) 


3C46 


0.437 


G 


<3.6 






3C268.1 


0.974 


G 


4.5 


1.2 


0.4 


3C268.3 


0.371 


G 


4.1 


0.8 


0.4 


3C268.4 


1.400 


Q 


6.1 


1.0 


0.6 


3C277 


0.414 


G 


2.5 


1.0 


0.2 


3C280.0 


0.998 


G 


12.0 


1.3 


1.2 


3C280.1 


1.659 


Q 


3.5 


1.4 


0.3 


3C286.0 


0.849 


Q 


428.0 


3.2 


40.0 


3C287.0 


1.055 


Q 


108.0 


1.4 


10.0 


3C288.1 


0.961 


Q 


4.4 


0.9 


0.4 


3C289.0 


0.967 


G 


3.4 


0.9 


0.3 


3C292.0 


0.713 


G 


3.4 


1.1 


0.3 


3C293.1 


0.709 


G 


< 3.0 






3C295.0 


0.461 


G 


35.6 


2.1 


4.0 


3C305.1 


1.132 


G 


3.4 


1.3 


0.3 


3C309.1 


0.905 


Q 


385.8 


2.5 


38.0 


3C313.0 


0.461 


G 


< 2.7 






3C320.0 


0.342 


G 


3.7 


0.8 


0.3 


3C322.0 


1.681 


G 


< 3.0 






3C323.0 


0.679 


G 


1.9 


0.8 


0.2 


3C325.0 


0.860 


G 


4.3 


1.0 


0.4 


3C337.0 


0.635 


G 


< 3.0 






3C343.0 


0.988 


Q 


6.6 


1.6 


0.6 


3C343.1 


0.750 


G 


8.9 


1.4 


0.8 


3C345.0 


0.593 


Q 


3480.0 


40.0 


300.0 








2150.0 


12.0 


100.0 


3C352.0 


0.806 


G 


5.8 


1.4 


0.6 


3C356.0 


1.079 


G 


< 3.0 







Uncertainties are given at la level 
upper limits at 3a 
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Table 2. IRAS Observations of 3C sources 



name 


12/im 


25 ptra 


60 /im 


100 iJ.rn 




(mjy) 


(mJy) 


(rnjy) 


(rnJy) 


3C46 


< 105 


< 66 


< 141 


< 350 


3C268.1 


< 60 


< 60 


< 75 


< 300 


3C268.3 


< 80 


< 74 


< 114 


< 345 


3C268.4 










3C277 


< 120 


< 100 


170±40 


< 450 


3C280.0 


< 66 


< 90 


90±30 


< 300 


3C280.1 


< 90 


< 75 


< 90 


< 390 


3C286.0 


< 75 


< 105 


< 90 


<270 


3C287.0 


< 100 


< 120 


< 150 


< 180 


3C288.1 


< 75 


< 60 


< 75 


< 180 


3C289.0 


< 75 


70±20 


100±40 


< 210 


3C292.0 


< 60 


< 75 


< 90 


< 240 


3C293.1 


< 120 


< 180 


< 143 


< 330 


3C295.0 


< 30 


< 60 


< 66 


< 210 


3C305.1 


< 60 


< 60 


< 60 


280±75 


3C309.1 


< 75 


80±14 


< 120 


< 600 


3C313.0 


< 70 


< 90 


170±32 


240±80 


3C320.0 


< 60 


< 66 


110±33 


270±90 




< DO 


< 4o 


90±2o 


zOO± 10 




< DD 


< 04 


< 9U 


< ioU 


3C325 

^w' <J fJ . W 


< 60 


< 45 


< 60 


< 180 


3C337.0 


< 60 


< 60 


< 80 


< 250 


3C343.0 


< 36 


< 45 


< 60 


< 180 


3C343.1 


< 45 


< 45 


< 75 


< 450 


3C345.0 


160±20 


310±20 


700±30 


1140±70 


3C352.0 


< 75 


< 75 


< 100 


< 240 


3C356.0 


< 75 


< 75 


< 75 


< 400 


averages RGs 


10±5 


< 15 


36±5 


59±30 


averages QSOs 


27.8±11 


23±8 


29±8 


62.6±27 



averages were computed with 



SUPERSCANPI procedures without 3C345 
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Table 3. ISO Observations of 3C sources 



name 


5/im 


7/im 


12/im 


25 




60 ^rn 


90 /im 


170 /im 


200 urn 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


3C46 












<50 


<40 


<300 


<480 


3C268.3 






<37 






<84 


<41 


<180 


<500 


3C268.4 














74±20 


430± 87 




3C280.0 


<7 


<25 


<90 


<; 


120 


120±40 


76±21 


<260 




3C286.0 


<10 










120±40 


100±30 






3C287.0 


<20 




<70 






<120 


<100 






3C288.1 






2.28±0.33'' 






<180 


<136 


<190 




3C295.0 


<20 


1.41±0.43'' 








160±40 


140±35 




<300 






















3C309.1 


20±8 


<30 


8.18±0.99<' 






100±30 


<180 


<260 




3C313 














<200 






3C325.0 


<50 


<20 


<60 






150±40 


100±30 


<100 




3C337.0 












<90 


<54 


<150 


<600 


3C343.0 






1.39±0.27" 














3C343.1 






0.80±0.27" 








<120 


220±40 


<250 


3C352.0 






<1.2 " 








<110 


135±70 




3C356.0 






0.83±0.28 " 








110±65 


<70 





" ISOCAM data 



